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The nonequilibrium dynamic phase transition, in the ki- 
netic Ising model in presence of an oscillating magnetic field, 
has been studied both by Monte Carlo simulation and by solv- 
ing numerically the mean field dynamic equation of motion 
for the average magnetisation. In both the cases, the Debye 
'relaxation' behaviour of the dynamic order parameter has 
been observed and the 'relaxation time' is found to diverge 
near the dynamic transition point. The Debye relaxation of 
the dynamic order parameter and the power law divergence 
of the relaxation time have been obtained from a very ap- 
proximate solution of the mean field dynamic equation. The 
temperature variation of appropiately defined 'specific-heat' is 
studied by Monte Carlo simulation near the transition point. 
The specific-heat has been observed to diverge near the dy- 
namic transition point. 

PACS number(s): 05.50.+q 



I. INTRODUCTION 

The dynamic response of the Ising system in presence 
of an oscillating magnetic field has been studied exten- 
sively [1-3] in the last few years. The dynamic hysteresis 
[1-3] and the noncquilibruim dynamic phase transition 
[4-6] are two important aspects of the dynamic response 
of the kinetic Ising model in presence of an oscillating 
magnetic field. This kind of phase transition in the Ising 
model was first studied by Tome and Oliviera [4]. They 
solved the mean field (MF) dynamic equation of mo- 
tion (for the average magnetisation) of the kinetic Ising 
model in presence of a sinusoidally oscillating magnetic 
field. They have defined the order parameter as the time 
averaged magnetisation over a full cycle of the oscillat- 
ing field and showed that depending upon the value of 
the field and the temperature, the order parameter takes 
nonzero value from a zero value. Precisely, in the field 
amplitude and temperature plane there exists a phase 
boundary separating dynamic ordered (nonzero value of 
order parameter) and disordered (order parameter van- 
ishes) phase. They [4] have also observed and located 
a tricritical point (TCP), (separating the nature (dis- 
continuous/continuous) of the transition) on the phase 
boundary line (see Fig.l). However, such a mean field 
transition is not truely dynamic in origin and exists even 



in the static (or zero frequency) limit. This is because, 
if the field amplitude is less than the coercive field (at 
temperature less than the transition temperature with- 
out any field) , then the response magnetisation varies pe- 
riodically but asymmetrically even in the zero frequency 
limit; the system remains locked to one well of the free 
energy and cannot go to the other one, in the absence of 
fluctuation. 

Lo and Pelcovits [5] first attempted to study the dy- 
namic nature of this phase transition in the kinetic Ising 
model by Monte Carlo (MC) simulation. Here, the tran- 
sition disappears in the zero frequency limit; due to 
the fluctuations, the magnetisation flips to the direc- 
tion of the magnetic field and the dynamic order pa- 
rameter (time averaged magnetisation) vanishes. How- 
ever, Lo and Pelcovits [5] have not reported any precise 
phase boundary. Acharyya and Chakrabarti [6] studied 
the nonequilibrium dynamic phase transition in the ki- 
netic Ising model in presence of oscillating magnetic field 
by extensive MC simulation. They [6] have also identi- 
fied that this dynamic phase transition is associated with 
the breaking of the symmetry of the dynamic hysteresis 
(to — h) loop. In the dynamically disordered (value of 
order parameter vanishes) phase the corresponding hys- 
teresis loop is symmetric, and loses its symmetry in the 
ordered phase (giving nonzero value of dynamic order 
parameter). They [6] also studied the temperature vari- 
ation of the ac susceptibility components near the dy- 
namic transition point. The major observation was that 
the imaginary (real) part of the ac susceptibility gives 
a peak (dip) near the dynamic transition point (where 
the dynamic order parameter vanishes). The conclusions 
were: (i) this is a distinct signal of phase transition and 
(ii) this is an indication of the thermodynamic nature of 
the phase transition. 

It is worth mentioning here that the statistical distri- 
bution of dynamic order parameter has been studied by 
Sides et al [7]. The nature of the distribution changes 
near the dynamic transition point. They have also ob- 
served [7] that the fluctuation of the hysteresis loop area 
becomes considerably large near the dynamic transition 
point. 

In the equilibrium critical phenomena, both the length 
and time scales diverge at criticality. This gives rise 
to the singularities in various thermodynamic quantities, 
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like the specific heat and the relaxation time. Can one ex- 
pect similar kind of features in the case of this noncquilib- 
rium dynamic phase transition problem ? To be specific, 

(i) Is there any 'relaxation time', in this nonequlibrium 
problem which diverges at the dynamic transition point ? 

(ii) Can there be any appropiately defined 'specific-heat', 
which will show singular behaviour at the transition point 
? 

The main motivation of this paper is to get some 
answers of these questions at least numerically. The 
nonequilibrium dynamic phase transition in the kinetic 
Ising model in presence of an oscillating magnetic field 
has been studied by MC simulation. Also the MF dy- 
namic equation has been solved numerically, to compare 
the results. The 'relaxation' behaviour (defined in the fol- 
lowing section) of the dynamic order parameter [4] and 
the behaviour of 'specific-heat' (defined in section III) 
near the dynamic transition point arc studied by MC 
simulation. It may be mentioned here that the prelim- 
inary results of 'specific- heat' singularity near the dy- 
namic transition point were reported briefly in ref [8]. 
More detailed results are reported here. The 'relaxation' 
behaviour has also been studied here by solving numer- 
ically the mean field (MF) dynamic equation of motion 
of the kinetic Ising model in presence of an oscillating 
magnetic field. The MF equation has also been solved 
exactly in the linearised limit and studied the 'relaxation' 
behaviour of the dynamic order parameter, near the dy- 
namic transition point. The paper is organised as fol- 
lows: In section II the 'relaxation' behaviour of the order 
parameter near the dynamic transition point has been 
studied both by Monte Carlo simulation and by solv- 
ing mean field dynamic equation of motion of the kinetic 
Ising model. In section III the temperature variation of 
the 'specific- heat' has been studied near the transition 
point only by Monte Carlo simulation. A brief summary 
of all the results is given in section IV. 

II. 'RELAXATION' BEHAVIOUR OF THE 
DYNAMIC ORDER PARAMETER 

A. Monte Carlo Study 

1. The Model and Simulation 

A ferromagnetically interacting (nearest neighbour) 
Ising model in presence of a time varying magnetic field 
can be represented by the Hamiltonian 

ff = -E ^fsj-M*)£s? (2-1) 

<ij> i 

Here, sf (= ±1) is Ising spin variable, Jjj is the interac- 
tion strength and h(t) = hocos(u)t) represents the oscil- 
lating magnetic field, where ho and u> are the amplitude 
and the frequency respectively of the oscillating field. 



The system is in contact with an isothermal heat bath at 
temperature T. For simplicity all Jij are taken equal to 
unity and periodic boundary condition is chosen. 

A square lattice of linear size L(= 100) has been con- 
sidered. At any finite temperature T and for a fixed 
frequency (oj) and amplitude (ho) of the field, the dy- 
namics of this system has been studied here by Monte 
Carlo simulation using Glauber single spin-flip dynamics 
with the Metropolis rate of spin flip. Each lattice site 
is updated here sequentially and one such full scan over 
the lattice is defined as the time unit (Monte Carlo step 
or MCS) here. The initial configuration has been chosen 
such that the all spins are directed upward. The instante- 
neous magnetisation (per site), m(t) = (l/L 2 ) Y\ sf has 
been calculated. From the instanteneous magnetisation, 
the dynamic order parameter Q = ^ § m(t)dt (time av- 
eraged magnetisation over a full cycle of the oscillating 
field) is calculated. 

2. Results 

Fig.l shows the schematic diagram of the dynamic 
phase boundary in the field amplitude (ho) and tempera- 
ture (T) plane. For small values of ho and T the dynamic 
order parameter Q is nonzero and the corresponding dy- 
namic hysteresis loop (m — h loop) is asymmetric and 
for larger values of ho and T the dynamic order parame- 
ter Q vanishes, correspondingly, the m — h loop becomes 
symmetric (inset of Fig.l). The dynamic transition tem- 
perature (Td) is a function of field amplitude (ho). The 
transition across the dotted line (in Fig.l) is discontinu- 
ous and that across the solid line is continuous. For very 
small values of ho the nature of the dynamic transition 
is continuous. In this paper, all studies are done in the 
region where Q undergoes always a continuous transition. 

It has been observed carefully that the dynamic order 
parameter Q does not acquire the stable value within 
the first cycle of the oscillating field. It takes several 
cycles (of the oscillating field) to get stabilised i.e., it 
shows 'relaxation' behaviour. Starting from the initial 
(all spins are up) configuration, the Q has been calculated 
for various number (say n-th) of cycles of the oscillating 
magnetic field and plotted (inset of Fig. 2) against the 
number of cycles (n). Each value of Q shown here has 
been obtained by averaging over 100 random Monte Carlo 
samples. Inset of Fig. 2 shows a typical 'relaxation' 
behaviour of the dynamic order parameter Q. This has 
been plotted for fixed values of to = 2n x 0.04, h = 1.0 
and T= 1.5. It shows that the dynamic order parameter 
Q is relaxing as the time (number of cycles) goes on. The 
best fit curves shows that the 'relaxation' is exponential 
type. So, one can write Q ~ Qoexp(— n/T), where T is 
the 'relaxation' time which provides the 'time scale' for 
this nonequilibrium problem. The physical interpretation 
of r is, the number of cycles required, so that Q becomes 
1/e times of its initial value (value of Q at starting cycle). 
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From the exponential fitting, the 'relaxation' time (T) has 
been meseared. The temperature (T) variation, for fixed 
values of oj and ho, of this 'relaxation' time T has been 
studied (in the disordered region of dynamic transition) 
and displayed in Fig. 2. The temperature (T) variation 
of r are shown (Fig. 2) for two different values of ho(= 
0.5 and 1.0) and for a fixed value of u> — 2ir x 0.04 here. 
From the figure (Fig. 2) it is clear that the relaxation time 
r diverges near the dynamic transition point (where Q 
vanishes) in the both cases (ho — 0.5 and 1.0). 



B. Mean field study 

1. Mean field equation of motion and numerical solution 

Although as mentioned earlier, the mean field system 
does not undergo a true dynamic transition (as the tran- 
sition exists even in the static limit), the mean field case 
has been considered here as a pathological one. 

The time evolution of the average magnetisation (un- 
der mean field approximation) in presence of an oscillat- 
ing magnetic field can be described by the equation 



dm 



h cos(ujt) 



dm 

t = —to + tanh 

dt 



m(t) + h(t)\ 



T 



(2.2) 



where m(t) is the instanteneous magnetisation, h(t) = 
hocos(ujt) is a sinusoidally oscillating magnetic field, T is 
the temperature and r is a constant. 

This equation has been solved by fourth order Runge- 
Kutta method taking r = 27rx0.01 and dt— 0.01. The ini- 
tial boundary condition is m(0) = 1.0. From the numer- 
ical solution for the instanteneous magnetisation m(t), 
the dynamic order parameter Q(= =r- § m(t)dt) has been 
calculated. 



dt T 
where e = 1 — 1/T. The solution of the above equation 
is 

m(t) = cxp(— ft It) + m cos(Lut — cf>), 

where mo and <f> are two constants. The value of the dy- 
namic order parameter Q at n— th cycle of the oscillating 
field is 

Q _ E. I m (t)dt = [ m(t)dt, 
2n J 2n Jt n -i 

where t n — 2im /u>. The value of Q, at the n-th cycle, 
can be written as 

Q = Qocxp( -) = Q exp(-n/F) 

TLO 

Qo is a constant independent of n. The above form shows 
that Q relaxes exponentially with the number of cycles 
(n) of the oscillating field. The 'relaxation' time T is 
equal to -^e -1 ■ It should be noted here that the dynamic 
transition occurs at T — 1 in the limit h — > [4]. So, 
for ho —> near the dynamic transition point (where the 
linearisation holds good) the behaviour of relaxation time 
is 

r-e- 1 ~ (T - T d (h - 0))- 1 

which shows the power law (exponent is unity) divergence 
of the 'relaxation' time at the dynamic transition point. 



III. BEHAVIOUR OF 'SPECIFIC-HEAT' NEAR 
THE TRANSITION POINT 

The time averaged (over a full cycle) cooperative en- 
ergy of the system may be defined as 



2. Results 

The inset of Fig. 3 shows a typical 'relaxation' of the dy- 
namic order parameter Q for uj = 2tt x 0.02, ho = 0.4 and 
T = 0.765. Here, also the exponential type of relaxation 
is observed and the 'relaxation' time has been measured 
in the same way, discussed earlier (in the MC case). Fig 3 
shows the temperature variation of the 'relaxation' time 
r for u — 2ir x 0.02 and two different values of ho (= 0.3 
and 0.4). Here also from the figure (Fig. 3) it is clear that 
the typical time scale or the 'relaxation' time T for this 
nonequilibrium problem, diverges, for both the cases (ho 
= 0.3 and 0.4), near the dynamic transition point (where 
Q vanishes). 



3. An approximate solution of MF equation 

In the limit of ho — > and T > 1, the equation ( |2.2| ) 
can be linearised (i.e., linearising tanh term) as 



(a) E c 



(uj/2ttL 2 ) I 



Z s i sZ j 
><*j> 



dt, 



and the time averaged (over a full cycle) total energy 
(including both cooperative and field part) of the system 
can be written as 



(b) E t , 



(uj/2-kL 2 ) 




dt. 



,<y> 



The temperature variations of E tot and E coop have been 
studied. The specific heat's are defined as Ctot — 
dEtot/dT and C coop — dE coop /dT. The temperature vari- 
ations of the 'specific- heat's have also been studied and 
found prominent divergent behaviour near the dynamic 
transition point (where Q vanishes). 

Here, again a square lattice of linear size L (=100) 
has been considered. Both E coop and E to t are calculated 
using MC simulation. Each data point has been obtained 
by averaging over 100 different MC samples. 
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A. Results 



The temperature derivatives of E coop and E to t can 
be denned as the 'specihc-heat's for this nonequilibrium 
problem. The temperature variations of Q , C coop (= 
dE coop / dT) and Ctot(= dE to t/dT) have been studied. 
The values of ho (=0.4 and 0.8) are chosen here in such 
a way that Q always undergoes a continuous transition. 
The temperature variations of Q, C coop has been shown 
in Fig. 4. Inset shows the variation of total coopera- 
tive energy E coop (per spin) with temperature (T). In 
this case the frequency (ui) of the field is kept fixed (w 
= 0.0628). Fig. 5 shows the temperature variation (for 
the same values of w, h and T) of C to t and the inset 
shows the temperature variation of the total (coopera- 
tive + field) energy (per spin). From the figure it is clear 
that the appropiately defined 'specific- heat' s C coop and 
Ctot diverge near the dynamic phase transition point. 



IV. SUMMARY 

The nonequlibrium dynamic phase transition, in the 
kinetic Ising model in presence of oscillating magnetic 
field, is studied both by Monte Carlo simulation and by 
solving the meanfield dynamic equation of motion. 

Acharyya and Chakrabarti [6] observed that the com- 
plex susceptibility components have peaks (or dips) at 
the dynamic transition point. Sides et al [7] observed that 
the fluctuation in the hysteresis loop area grows (seems 
to diverge) near the dynamic transition point. 

In this study it is observed that the 'relaxation time' 
and the appropiately defined 'specific-heat' diverge near 
the dynamic transition point. All the results are obtained 
here numerically. No attempts were made to extract any 
exponent values from the numerical studies. 

It should be mentioned that recent experiments [9] on 
ultrathin ferromagnetic Fe/Au(001) films have been per- 
formed to investigate the frequency dependence of hys- 
teresis loop areas. Recently, attempts have been made 
[10] to measure the dynamic order parameter Q experi- 
mentally, in the same material, by extending their previ- 
ous study [9]. The dynamic phase transitions has been 
studied from the observed variation of Q. However, the 
detailed study of the dynamic phase transitions by mea- 
suring variations of associated response functions (like 
the ac susceptibility, 'specific-heat', correlations, relax- 
ations etc) have not been studied experimentally. 
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Figure Captions 

Fig. 1 Schematic diagram of the dynamic phase boundary 
in the field amplitude (ho) and temperature (T) plane. 
The dotted line is the boundary of the discontinuous 
transition and the solid line represents the boundary of 
continuous transition. The small circle represents the tri- 
critical point (TCP). Insets demonstrate the breaking of 
the symmetry of the dynamic hysteresis (m — h) loop due 
to dynamic transition. 

Fig. 2 Monte Carlo results of the temperature (T) vari- 
ation of 'relaxation' time (r) for two different values of 
field amplitudes (ho):the bullet represents ho = 1.0 and 
the diamond represents ho — 0.5. Solid lines show the 
temperature (T) variations of dynamic order parameter 
Q. Inset shows a typical 'relaxation' of Q plotted against 
the number of cycles (n). The solid line is the best fit ex- 
ponential form of the data obtained from MC simulation. 
Here, L = 100, uj = 2tt x 0.04. 

Fig. 3 Mean field results of the temperature (T) variation 
of 'relaxation' time (r) for two different values of field 
amplitudes (ho): the filled triangle represents ho = 0.4 
and filled square represents ho = 0.3. Solid lines represent 
the temperature variation of the dynamic order param- 
eter Q. Inset shows a typical 'relaxation'of Q plotted 
against the number of cycles n. The solid line is the best 
fit exponential for of the data obtained from the solution 
of equation 2.2. Here, ui = 2ir x 0.02. 



Fig. 4 Monte Carlo results of the temperature variations 
of C coop for two different values of ho: the filled square 
represents h = 0.8 and the filled triangle represents ho 
= 0.4. Solid lines represent the temperature variations 
of Q. Inset shows the temperature variations of E coop for 
two different values of ho'.il) ho — 0.8 and (II) ho = 0.4. 
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Here, L = 100, uj = 2tt x 0.01. 

Fig. 5 Monte Carlo results of the temperature variations 
of Ctot for two different values of ho'- the filled square 
represents h = 0.8 and the filled triangle represents ho 
= 0.4. Solid lines represent the temperature variations 
of Q. Inset shows the temperature variations of E to t for 
two different values of ho: (I) ho = 0.8 and (II) ho = 0.4. 
Here, L = 100, u = 2tt x 0.01. 
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